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We extract the P-wave %% phase shift for five values of pion relative momenta, which gives infor- 
mation on the p resonance. The Breit-Wigner formula describes the n% phase shift dependence 
nicely and we extract m p = 792(7)(8) MeV and the coupling gp nK — 5.13(20) at our m n = 266 
MeV. We extract the P-wave scattering length af^ = 0.082(10)(3) fm 3 from the state with the 
lowest pion relative momenta. 

We also determine the S-wave p% phase shift for two values of relative momenta, which provides 
parameters of the lowest axial resonance a\ (1260). Using the Breit-Wigner fit we extract m a \ = 
1 .44(4) GeV and the coupling g U[ pn — 1 • 1 (3) GeV. From the lowest state we also extract the pn 
scattering length af= = 0.23(12) fm for our m n . 

The simulation is performed using one Nf = 2 ensemble of gauge configurations with clover- 
improved Wilson quarks. The phase shifts are determined from the lowest two energy-levels, 
which are obtained by the variational analysis with a number of quark-antiquark and meson- 
meson interpolators. The correlation functions are calculated using the distillation method with 
the Laplacian Heaviside (LapH) smearing of quarks. 
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1. Introduction 

Extracting the width of a hadronic resonance R from lattice QCD is challenging. The only 
proper method used up to now applies to resonances R that appear in the elastic scattering of 
two hadrons H1H2 H\H2- First the elastic phase shift 8{s) for //1//2 scattering has to be 

determined from the lattice for several values of s = E^ M = E — P , where E and P are the energy 
and the total momentum of the H1H2 system. Liischer has shown that the energy E of two hadrons 
in a box of size L ~ few fm provides the value of the infinite- volume elastic phase shift 8 (s) at 
s = E 2 — P 2 [p. His relation between 5 and E for P = was generalized to P / in [g |, §]. In 
practice, one or two lowest energy levels E are extracted and a few choices of P are used in order 
to extract 5 (s) at different values of s = E 2 — P 2 . 

The resulting 8 (s) can be fit with a Breit-Wigner (or any other desired) form, where both are 
related via the scattering amplitude a\ for the /-th partial wave 

-V~sr R (s) e 2i5 ^ - 1 2 2 p* 2l+l 

«/ = 2 , ■ rr < \ = ^ ° r Vsr R {s) cot S{s) =m R -s, r R (s) °c g RHH 

s — m R + i^JsT R (s) 2i s 

(1.1) 

and /?* is the momentum of Hi and H2 in their center-of-momentum (CMF) frame. This relation 
can be used to extract the mass m R and the width F R = T R {m\) of the resonance from lattice data 
on 8{s). The width depends significantly on the phase space and therefore on m n , so it is common 
to extract the coupling gRHtfy, which is expected to depend only mildly on m n . 

Among all the meson resonances only the p meson width has been determined properly using 
this method. The first lattice determination was done by PACS-CS in 2007 [j|]. Since then, several 
studies of the p have been carried out [^, ^J, with the most recent ones [||, ||, [T(J. In this talk we 
present our recent study of the p |9||, which achieves the smallest statistical errors (on one ensemble 
only, however) on the resulting 8(s), m p and T p due to several improvements listed below. 

We also extract the S-wave pn elastic phase shift, which enables us to extract the mass m a \ 
and the width r a i of the lowest lying axial resonance a\ (1260). The lattice study of this resonance 
is especially welcome as the experimental knowledge on it is very poor: the width has a wide range 
T^f = 250 — 600 MeV [jllj], and none of its branching ratios have been reliably determined 1 [jnj]. 
To our knowledge, this is the first lattice study aimed at the pn scattering and F a \. 



2. Lattice simulation 

We use 280 Nf = 2 configurations with tree-level clover-improved Wilson dynamical and 
valence quarks, corresponding to m n a = 0.1673(16) or m n = 266(3)(3) MeV [jl2|]. The lat- 
tice spacing a = 0.1239(13) fm was determined using the Sommer parameter ro []|] and our 
N[ x Nj = 16 3 x 32 is rather small, allowing us to use the powerful but costly full distillation 



method [13]. We combine periodic and anti -periodic propagators in time to reduce the finite Nj 
effects [Hi. 



All final states are quoted just as "seen" in [11 IP 
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Figure 1: Contractions for our p and a\ correlators with qq and meson — meson interpolators (7 = 1). 

3. p resonance and tztz phase shift 

The details of our lattice simulation aimed at %% phase shifts and the p resonance have been 
published in [0]. In this talk, we emphasize the most important steps and results. 



The n + n 



n + n scattering is elastic below the 4n threshold yfs < 4m n and we can 



apply Liischer's method. We determine the lowest two energy-levels of the p n + 7i coupled 
system with J PC = 1 and = | 1,0) for the following cases of total momentum P 
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and all permutations in direction P and p -polarization. We display the symmetry group, the irre- 
ducible representation and the decay mode, which applies to three cases of P ^ [8], |9|. 

Other simulations aimed at F p used at most one quark-antiquark interpolator and one %% 
interpolator for each P. We use 15 quark-antiquark interpolators &*~^"^ w and one nit interpolator 
for each P, where each pion is projected to a definite momentum: 



'i=l...,5 



£ ^ [u s (x) &i e iPx u s {x) - d s {x) &i e iPx d s (x)\ {s = n,m,w) , (3.1) 

X 

m) = ^ + (pi)7T ( P2 ) - 7r-(pi)7r + (p 2 )] , ^(Pi) = ^n{x)y 5 ^e^q n {x) . 

x 

Quark-antiquark interpolators have five different color-spin-space structures The quarks are 



smeared using the Laplacian Heaviside (LapH) smearing proposed in []13j], i.e., 
^©(a, 2 + V 2 )4=£0(a, 2 + A«; 



v (*) v (*)t. 



s = n {narrow), m {middle), w {wide) , (3.2) 



k=i 



where different truncations N v = 96, 64, 32 correspond to three different widths s = n,m,w |]9p. 

The 16 x 16 correlation matrix Cij{tf,t{) = (0| $i{tf)v:{ti)\0) necessitates the inclusion of the 
contractions in Fig. [j] The contractions were computed using the full distillation method, which is 



based on the LapH smeared quarks ( ft.2| ) []13J and leads to relatively precise results for all types of 
contractions in Fig. [jj We apply this method for the first time to extract a meson width. We also 
propose how to apply it for interpolators with different smearing widths in the same variational 
basis [gj. All correlators are expressed in terms of the so-called perambulators in Appendix A of 
[g]. The resulting correlators are averaged over all source time-slices over all directions of P and 
p polarization. 

The lowest two energies of the system are determined using the Generalized Eigenvalue Method 



(GEVP) [14] and the dependence on the choice of the interpolators in the variational basis is ex- 
plored in [B]. The lowest energy level is robust to this choice. We find that the first excited energy 
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to 




s s 

Figure 2: The phase shift 8 (in degrees) for %% scattering in P-wave and ((op*) 3 /Vsa 2 ) cot 8 as a function 
of s, together with a Breit-Wigner fit. 



level cannot be reliably obtained without the kk interpolator in the basis, and that more than two 
interpolators are required at least in the case P = ^(1,1,0). The extracted six energy levels for 
our preferred interpolator choice [[J are given in Table III of 

Each of the six energy levels gives the value of the phase shift 8 (s) at s = E 2 — P 2 ( 2 ) via the 
Liischer formula for P = or its generalization to P 7^ [g, |J]. We independently confirmed the 
needed relations and compiled them in One of these levels, i?2(P = 0), is above the inelastic 
threshold y/s > 4m n and we omit it from further analysis. 

The resulting phase shifts for five different values of s are plotted in Fig. |2[ The phase 
shift has relatively small errors and exhibits a resonating behavior, which allows us to extract 
nip and T p or rather the coupling g pnK . We use the Breit-Wigner relation ( |0| ) together with 
r ( J ) = gpTtK P* 3 /{6xs), which leads to 



V* Son a 



L— cot 8 (s) = -p- (m 2 p -s) . (3.3) 



This allows a linear fit in s (Fig. |||) to extract m p and g p7lK given in Table p]. The resulting m p a = 
0.4972(42) is slightly lower than the naive value m p mve a = 0.5107(40), which is extracted from 
the ground state with P = 0. We also extract the P-wave scattering length af^ { = 0.082(10) (3) fm 3 



(defined as a, = lim s ^ 8(p*)/p* + [(Bp) from the state with the lowest 3 p*a = 0.1076(36) and 
8 = 3.03(6)°. This qunatity is not directly measured, so we compare it to the typical value af^ ~ 
0.038(2) (m^ /l> )~ 3 obtained by combining experiment and ChPT or Roy equations [[IB]]. 

A comparison of the resulting m p and g p7[7[ to two recent lattice simulations [|| 10] is com- 
piled in JlC[] - The Nf = 2 simulation with twisted mass quarks [^] and the Nf = 2+1 simulation 
with Wilson quarks were done at four/two values of m n and explicitly demonstrate the mild depen- 
dence of g p7[K on m n . All three results on g pn]t are relatively close to each other and close to the 
g C p7C7i = 5-97 extracted from Tp* 11 ' . The resonance mass m p of ^ is ~ 11% higher than ours, while 
m p °f [§] is — 21% higher than ours, at comparable m n . Note that all three simulations get the 
resonance m p within 3% from the value of m p mve , which implies that the simulations differ already 
in m' p aive . Possible causes for different m' p " ve could be discretization effects or scale fixing of all 



2 We use the discrete dispersion relation cosh(^fsa) = cosh(£a) — 2£^ =1 sin 2 (/ > j : a/2) instead of the continuum one 
= E 1 — P 2 to analyze the p g], We analyze the a\ using the continuum dispersion relation. 



3 The next state leads to affj consistent with the value obtained from the lowest state. 
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m p [MeV] g pn% 


n 7171 
a l=i 


m al [GeV] 


ga x pn [GeV] 


«f=o [ fm l 


latt 


792(7)(8) 5.13(20) 


0.082(10)(3) 


1.44(4) 
1.43(5) 


1.1(3) 
1.7(4) 


0.23(12) using mp 
0.56(23) using m n p aive 


exp 


775.5 5.97 


0.108(5) * 


1.23(4) 


< 1.35(30) 


not meas. 



Table 1: Our lattice results for the resonance properties [[J, compared to the experimental values. The 
results related to a\ depend on the choice of the input p mass: m p or m' p " ve . The experimental value of af^ 
is obtained combining experiment with ChPT or Roy equations. 



three simulations, flavor breaking of twisted mass quarks or partial quenching of the strange 
quark [Q, H]. Additional causes for the different m p values could be the small interpolator basis in 



10[ ] or tne small box L ~ 2 fm in [9fl. The exponentially suppressed terms, which are neglected 



in Liischer formulae, may not be completely negligible for our L ~ 2 fm, which is a systematic 
uncertainty of our simulation. We are planing a simulation at larger L to explore possible finite size 
effects. We believe, however, that our small L does not influence our m n p mve , as the first excited state 
n(2n/L)n(-2n/L) at P = hardly affects the m n p aive ground state. 

Our 8 (s) agrees reasonably well with the prediction of the lowest 4 order of Unitarized Chiral 
Perturbation Theory [ff™, which has been recalculated for our m n = 266 MeV 



4. The pn phase shift and a\ resonance 

We study the S-wave scattering of pn, where the resonance «i (1260) appears, for the total 
momentum P = 0. The scattering is elastic at least until a\ (1260) on our lattice since K*K cannot 
be created on our Nf = 2 ensemble. The ground scattering state is p (O)tt(O) in the non-interacting 
limit. The scattering particle p (0) is almost stable on our lattice, since its lowest decay channel 
7i(2n/L)7i(—27i/L) is significantly higher in energy. 

We use 9 quark-antiquark interpolators ^~ 1 ' v "' w and one p (0)7r(0) interpolator, all with J PC = 
1++, 1^/3) = 1 1,0) andP = 0: 

= L v5 A M >* e iPx u s (x) - {u s <-> d s } (s = n,m,w) , (4.1) 

x,i 

= L 73 "*(•*) w Ai V ^ eiPx m' m *( x ) - {»* ° d *} ( s = n ^ m M » 

^(0 = £ 75 "s-W AiYj ^[e iPx V, - V/e iP > 4 .(.x) - {u s o rf,} (j = n,m,w) , 

^(0 = j 1 [n + (Q)p-(0)-n-(0)p + (0)} , ^(0) = £^5^ , p ± (0) = ^^A^-t^ , 

X X 

where V is the covariant derivative. The contractions in Fig. [I] are calculated using the full distilla- 
tion method and averaged over all source time slices ti and all a\ polarizations A. 

The effective mass for the lowest two eigenvalues are shown in Fig. |3| and the resulting E and 
p* are given in Table ^[ The CMF momentum p* is extracted using E = \J p* 2 + m| + \J P* 2 + m 2 p : 

4 One cannot make a fair comparison between our lattice result and the NLO prediction, since it depends on a number 
of LECs, and some of them have been fixed using m p from another lattice study, which gets a significantly higher m p . 
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Figure 3: The effective mass for lowest two eigenvalues in a\ channel (left). The combination p* cot 8/yfs 
as a function of s = E 2 — P 2 , where 8 is p 7T phase shift in S-wave (right). 



it is imaginary for the ground state below m n + m p threshold, and real for the first excited state. We 
take two choices for the input p mass: our main results are based on the resonance mass m p (green 
lines in Fig. ||), while m n p mve is taken for comparison (blue lines in Fig. ||). 

The S-wave phase shift 8 for P = is extracted using the well known Liischer relation ^ 

2 _ /. , „*t ,t\ p*^o 1 



p* cot 5 = 



pK 

2 l=0 



(4.2) 



which applies above and below threshold. The results are compiled in Table || for both choices of 
p mass. The first excited level gives 8 r* 90°, so it is sitting close to the top of the a\ resonance and 
m a \ w E2 holds. The ground state with imaginary p* gives imaginary 8, but the product p* cot 8 is 
real since Zqo(1, (^) 2 ) is reai - 

We parametrize r a i (s) = g^ p7[ p* /s and apply the Breit-Wigner relation ( |1.1[ ) to get 

1 



^cot8(s) 



(4.3) 



which applies in the vicinity of the resonance above or below threshold. Given the values of p* cot 8 
at two different values of s, we apply a linear fit (4.3) in s (shown in Fig. |3|) to extract m a \ and 
gcupn- The results are compiled in Table [j]. Our m a \ at m n = 266 MeV is about 14% higher than the 
experimental resonance a\ (1260). The first lattice result for g ax pn is valuable, since this coupling is 
not known experimentally. None of the a\ branching ratios have been measured, so we provide only 
the upper limit for g^ p % resulting from the total width F^f = 250 — 600 MeV. Our lattice result 
gaipn = 1-1(3) GeV is in agreement with the value g P a x e p n % ~ 0.9 GeV obtained using Unitarized 
Effective Field Theory approach Jl7| ] and converted to our convention. We extract also a^ Q from 
the ground state, which is sufficiently close to the threshold. The scattering experiment cannot 
be carried out since p is a quickly decaying particle, so we compare our a^ (m 7[ = 266 MeV) = 
0.23(12) fmtotf* (m p J' y 



0.37 fm obtained from Unitarized Effective Field Theory [18] 



5. Conclusions 



The lattice extraction of the phase shifts for elastic scattering has recently become possible 
also for the attractive resonant channels. We simulated the scattering in the p and a\ channels and 
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level fit Ea = yfsa p*a 8 p*cos(8)/y/s 



1 7-10 0.6468(73) 


i 0.065(13) 


i 7.1(54)° 


0.82(44) 


(using ntp) 




i 0.086(9) 


i 23(14)° 


0.34(14) 


(using m n p aive ) 


2 6-9 0.897(13) 


0.280(10) 


83.7(59)° 


0.034(33) 


(using ntp) 




0.272(10) 


88.9(59)° 


0.005(31) 


(using m n p aive ) 



Table 2: The results for the a\ <H> %p coupled channel with interpolators ff" 2 4 and GEVP reference time 
?0 = 5. The ground state is below pit threshold, so p* and 8 are imaginary, while the p* cot 8 is real. 

extracted the mass and the width of these two resonances as well as the scattering lengths in the 
corresponding meson-meson channels. 
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